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Abstract 

We reconsider the mass spectrum of double sine-Gordon theory where recent 
semiclassical results called into question the previously accepted picture. We use 
the Truncated Conformal Space Approach (TCSA) to investigate the claims. We 
demonstrate that the numerics supports the original results, and strongly disagrees 
with those obtained from semiclassical soliton form factor techniques. Besides the 
numerical analysis, we also discuss the underlying theoretical arguments. 

1 Introduction 

Double sine-Gordon theory has attracted interest recently chiefly because it is a prototype 
of non-integrable field theory which can be understood by application of techniques devel¬ 
oped in the context of integrable field theories [1]. It has several interesting applications 
such as to the study of massive Schwinger model (two-dimensional quantum electrodynam¬ 
ics) and a generalized Ashkin-Teller model (a quantum spin system) which are discussed 
in [1], Another application to the one-dimensional Hubbard model is examined in [2] 
(together with the generalized Ashkin-Teller model mentioned above). A further poten¬ 
tially interesting application of the two-(and multi-)frequency sine-Gordon model is for 
ultra-short optical pulses propagating in resonant degenerate medium [3]. 

In an earlier work the authors (and coworkers) studied this model using non-perturbative 
finite size techniques and established the phase diagram of the model [4] (which was recently 
extended to the multi-frequency generalization [5]). In particular we used form factor per¬ 
turbation theory introduced by Delhno and Mussardo in [1] to predict mass spectra and 
verify them using the truncated conformal space approach (TCSA) [6]. 
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Recently Mussardo et al. [7] applied a semiclassical soliton form factor technique de¬ 
veloped by Goldstone and Jackiw [8] to study the mass spectrum of the theory. They 
obtained results that contradict explicitely some of the results obtained using form factor 
perturbation theory in [4]. In principle, the truncated conformal space approach applied 
in [4] could decide the question, but the accuracy achieved there is not sufficient for this 
purpose. 

The present work reports a systematic study of this issue using improved TCSA nu¬ 
merics which do have the required accuracy. Besides the numerical comparisons we also 
address the theoretical arguments raised in [7] to support their case. We find a convincing 
and consistent picture, both numerically and theoretically, which allows us to clarify the 
issue. 

We would like to stress that the problem at hand is much more than a simple nu¬ 
merical discrepancy: it is the validity of the form factor perturbation theory and that of 
semiclassical techniques that is at stake here, both of which are very promising tools for 
non-pertnrbative investigation of non-integrable quantum held theories. 

We start with recalling some necessary facts about the double sine-Gordon model to 
set the stage in Section 2. Section 3 describes the mass spectrum obtained from form 
factor perturbation theory in [4] on the one hand and from the semiclassical soliton form 
factor technique on the other [7]. Section 4 describes the numerical work and the results 
obtained, while Section 5 is devoted to the underlying theoretical issues. We conclude 
briefly in Section 6. 


2 The double sine-Gordon model 


The classical model is defined by the Lagrangian density 


-^classical = -d^ipd^tp ~ /J 0 COS (3(fi - A 0 COS 



For Ao = 0 we recover standard sine-Gordon theory, where the spectrum is known to consist 
of a soliton doublet of mass 


M 0 = 



P 


and a continuum of breathers 


7T 

m = 2 Mn sin —a 


0 < a < 1 


We define the quantum version of the model using perturbed conformal held theory as in 
[4], The Hamiltonian can be written 

77pcft = H c= i + // j dx : cos flip : +A J dx : cos ^ : 
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where 


H c= i = j dx + \ (d x (p ) 2 ■ 

is the Hamiltonian of a c = 1 free boson CFT. The quantum spectrum for A = 0 consists 
of a soliton doublet of mass M which is related to the coupling p by [9] 

2r(A) A p 

" rf(l-A)^ 2r(^) ) ' Srr 

which satisfies the semiclassical asymptotics 

/i ~ ^M 2 , (3 ^ 0. 

which means that p can be identified with /i 0 in that limit (similar results hold for A and 
Ao). In the quantum theory, the breather spectrum is discrete and their numbers depend 
on the frequency (3 : 


m n 


7T n/3 2 

2M sin — - 


2 8t r-f3 2 


n = !,•••, 


8it 

w 


-1 


( 2 . 1 ) 


Double sine-Gordon model can be considered as a non-integrable perturbation of the 
integrable sine-Gordon held theory [1], Because of the presence of the second cosine term, 
the held must be taken to have the periodicity 


47T 



which is called a 2-folded model following [10]. Therefore at A = 0 the theory has two 
degenerate vacua corresponding to the classical solutions 

p — 0 and 


As a result, each breather comes in two copies, and soliton states must be labeled with the 
vacua, between which they mediate, thereby implementing a so-called ’kink’ structure [10]. 
In this paper we treat A as the non-integrable perturbing coupling constant, although it is 
equally possible to swap the roles of the cosine terms in this respect. 

The most interesting case is 8 — | when for small A the two vacua remain degenerate 
until a critical value is reached where the theory undergoes a phase transition. This 
was shown at the classical level (mean held theory) by [1] to be of second order in the 
Ising universality class, and it was established that this holds even after taking quantum 
fluctuations into account (non-perturbatively) in [4], This is the model in which we are 
the most interested, although we shall also use the 5 = 0 case as a testing ground for our 
numerical calculations. 
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3 Mass spectrum predictions 


3.1 Form factor perturbation theory 

Form factor perturbation theory (FFPT) was introduced in [11] to provide a way of ex¬ 
panding physical quantities in the non-integrable perturbing coupling. It was applied to 
the double sine-Gordon model in [4]. For our present purposes, we only need the mass 
shift of the first breather, which reads 1 


Sm m = ,_i_ r*> dt 

M l 71 Jo 


r i 

rirp 

t \ 

— 

/ dt —— > 

l 7r 

Jo 

sm t) 


where k — 1,2 labels the two copies of the first breather, p = (3 2 / (87t — /? 2 ), 

X/5FT (s^)' 


6M = 


x exp 


2 r (-££-) 

ZL \Stt-I 32 ) 

r°° dt 

lo T 


si 

47r 


sinh 2 | 

CtfO 

1 

2 sinh | 

im 

1 sinh (t) cosh ( 

h-S) 

u) 


and 


M = exp < 4 


dt sinh(t) sinh (pi) sinh((l + p)t) 




-2 t 


t 


sinli 2 1 


It is an important prediction of FFPT that the mass corrections vanish when 5 = 
that also 

5m^ = —Sm^ 


in general. 


(3.1) 


(3.2) 


(3.3) 
| and 


3.2 Masses from semiclassical soliton form factors 

In [7] a semiclassical method, based on the Goldstone- Jackiw semiclassical kink form factor 
[8], was used to determine the spectrum of double sine-Gordon theory. Here we give a 
brief summary of their results for the case 5 = |. They claim that the perturbation by 

cos f | j induces a linear mass correction to the mass of the nth breather, and obtain 

the following result: 

m "“ 1=m ” * 2 ptI si “ (4) ■ is cos (4 )) + ° G) (3 - 4) 

where L and S index the splitting up of the nth breather level into what they call a ’short’ 
and a ’long’ breather. 

1 We corrected some misprints in the published version of [4], 
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The main problem and the central issue in this paper, is that this disagrees with the 
FFPT results in the previous subsection that give a vanishing linear correction to the 
breather masses, and it also disagrees with the numerical studies in [4] where the masses 
of the two copies of any breather levels remained degenerate as long as A was lower than 
the critical value, where a second order phase transition induced a mass gap between the 
ground states and also lifted the degeneracy of the breather masses. However, the numerical 
studies in [4] were not aimed at studying the breather levels and the data produced for 
that work, while still accessible, have insufficient precision to decide whether a splitting 
of the order of (3.4) is actually present. In addition, the authors of [7] presented some 
arguments to support the result (3.4), to which we shall return later. For the moment the 
issue is to decide which of the predictions is borne out by the dynamics of the model. 

Using the soliton mass M 0 = 8y/JIo//3, the mass correction can be written as 


Sm n' S> -MR V / 1 . ( P 

M 0 M 2 \p 2 \ 


2 ' 

16 


n (3 
— cos ri- 

lG \ 


2 ' 

16 


7rrr/3 4 Ao 
768 Ml 


+ 0 


In terms of the dimensionless coupling constant 

-^o 


Xo = 


MS 


this can be written 


8 m ^ ,s ' > 7m 3 /? 4 




M 0 ~ 768 /vu ' ” V~ J (3 ' 5) 

where we again absorbed the difference between the classical and quantum mass/coupling 
relations into the higher order correction term. 

The dimensionless coupling of the quantum theory is defined by 


X = 


A 


9 P 
]\/[ 16 - 71-2 


(where the difference in the exponent of M is just the anomalous dimension of the per¬ 
turbing operator). In the semiclassical limit x only differs from xo by terms higher order 
in (3 2 as discussed in the previous section, and similarly for M and My. Therefore we can 
still use formula (3.5) to leading order in (3 in the quantum theory just replacing Mo by M 
and xo by x 

Sm^ _ to 3 /? 4 

M 768 X 

For our purposes we are only interested in the first breather (n = 1) and the comparisons 
to numerical results will be done at the following couplings: 



YT '■ ^m~ = ±0 - 0165367 * + • • • 

IT : ~M~ = ±0 - 0275T51 X + ■■■ (3.6) 
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4 TCSA analysis 

We now aim to get the mass spectrum solving the model directly using the so-called 
truncated conformal space approach (TCSA) introduced by Yurov and Zamolodchikov in 
[6]. It was extended to perturbations of c = 1 free boson models in [12] and applied to the 
double sine-Gordon model in [4]. 

The application of TCSA to perturbations of c = 1 free boson CFT is described in detail 
in [4] and references therein, therefore we do not go into the details here. We mention only 
that this method gives the finite volume spectrum of the theory, from which the masses can 
be extracted by locating the so-called scaling region in which the mass gap function (the 
energy difference between the one-particle state and the corresponding vacuum state) is 
approximately constant. More precisely, we locate the volume range where the magnitude 
of the numerical derivative of the mass gap function is smallest and define the mass as the 
value the mass gap function takes there. The Hilbert space can be split up into sectors 
differing in the value of the winding number of the held ip and the eigenvalue of total 
spatial momentum, and the mass gap functions we are interested in can be measured from 
the sector with no winding and zero momentum. 

We improved on the precision of the TCSA significantly with respect to [4], First, we 
chose smaller values for f3 which makes TCSA converge faster in the models at hand. This 
is also advantageous since it means we are closer to the classical limit and so we expect 
to be able to test the semiclassical prediction (3.4). Second, we pushed the cutoff much 
higher using more machine power, and also by applying the parity symmetries to project 
the Hilbert space into even and odd subspaces, using that the TCSA Hamiltonian is block 
diagonal in the appropriate basis. This made possible treating more than 18000 states in 
reasonable computing time, and allowed us to achieve the accuracy which makes it possible 
to judge the correctness of the semiclassical prediction (3.4). 

4.1 5 = 0: testing the numerics 

In this case the model has a symmetry 

<P “»■ (4- 1 ) 

which can be used to project the TCSA Hilbert space onto an even and an odd sector. 
This reduces the number of states and helps to achieve higher values of the cutoff. In 
order to be on safe ground, the program was run both with and without a projection for 
several values of /3, x an d the cutoff. The resulting spectra were always identical (within 
the precision of floating point arithmetics), independently of the values of the parameters 
and of the cutoff, which confirms that the symmetry projection was correctly implemented. 

We also compared the y = 0 to predictions of exact S matrix theory, which they 
matched perfectly within the truncation errors. The vacua can be found in the even, while 
the one-particle states in the odd sector. They can be matched against each other using the 
fact that the bulk energy constant of the two vacua are different, therefore the two vacua 
(and also the two corresponding one-particle states) have a different slope in the scaling 
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region. The redefinition A —> —A leaves the total spectrum invariant (in both the odd and 
the even sector separately), but exchanges the role of the two sets of states, therefore it is 
expected that 

mj 1,3 ' (—A) = mf' 1) (A) 

We choose the following value for the coupling 


P = 



1.4179631... 


For mass determination, the cutoff was placed at level 18 resulting in 9334 states in the 
even and 9319 states in the odd sector respectively. 

At this value of [3 there are 


breathers in the A = 0 spectrum which shows that this model is deeply in the semiclassical 
regime. The FFPT formula (3.1) gives the following prediction for the mass corrections: 


Sm^’ 2 ^ 

M 


±0.95392713... x x 


(4.2) 


We summarize the measured values of the masses in table 4.1 and plotted them against 
predictions of form factor perturbation theory in figure 4.1. It is obvious that the measured 
masses follow well the predictions from form factor perturbation theory, but there are 
deviations that grow with the coupling. These are expected to be higher order effects. 
However, second order corrections should cancel from the difference, since sending A to 
—A exchanges the identity of the two particles [4], In fact, the mass differences agree with 
form factor perturbation theory within the estimated TCSA errors which can be estimated 
from cut-off dependence of the energy levels and the flatness of the gap function in the 
scaling region. We can also fit the coefficient of (4.2) against the mass splitting measured 
in TCSA, for which we find 

0.95441 ± 0.00032 


The central value differs by only 0.0005 from (4.2), consistent with the la error range of 
the fit. 

It is also interesting to examine the mass correction in the classical limit. Calculating 
the curvature of the field potential around its minima at ip = 0 and ip = ^ gives 


Sm ( i ,2) _ A 0 

M 0 ~ Xo ’ X ° “ Mq 


(4.3) 


We can see that quantum corrections to the coefficient are of the order 0.05, and the TCSA 
method is able to measure them with one percent accuracy. It is also useful to remark that 
the semiclassical expansion of quantum corrections to mass ratios goes in the parameter 


87T 
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(a) The two masses, denoted here by ml and m2, from TCSA and FFPT 
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(b) The mass splitting from TCSA and from FFPT 

Figure 4.1: Comparing masses measured from TCSA and predicted using FFPT for f3 = 
,5 = 0. All data are in units of the unperturbed soliton mass M. ml and m2 denote 
the two masses m[ l '~\ respectively. The error bars of TCSA data are too small to visualize. 
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0.001 

0.27325 
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0.00381 

0.00380 

0.003 
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0.26947 

0.00572 

0.00572 

0.004 

0.27610 

0.26847 
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0.26851 

0.00763 

0.00764 

0.005 

0.27702 

0.26749 

0.27710 

0.26756 

0.00953 

0.00954 

0.006 

0.27796 

0.26652 

0.27805 

0.26661 

0.01144 

0.01144 

0.007 

0.27888 

0.26553 

0.27901 

0.26565 

0.01335 

0.01336 

0.008 

0.27982 

0.26455 

0.27996 

0.26470 

0.01527 

0.01526 


Table 4.1: Comparing masses measured from TCSA and predicted using FFPT for (5 = 
,5 = 0. All data are in units of the unperturbed soliton mass M. TCSA mass 
measurements have an estimated error of order ±0.00002. 


(as can be seen e.g. by expanding the unperturbed sine-Gordon breather/soliton mass 
ratio (2.1) for small 8) and for this particular value of 8 we have 


8tt 


0.08 


which explains why the corrections to the classical limit are so small. 


4.2 6 = 7r/2: the acid test 

In this case there is a symmetry 

2vr 

^ ( 4 - 4 ) 

which we used for projection. We again performed runs with and without a projection to 
confirm that the program was functioning correctly. In this case, both the even and the 
odd sector contains a vacuum and a one-particle state. 

For the value of 8 used in the previous section, we obtain the masses listed in table (4.2). 
At the same truncation level (18) the dimension of the even sector is 9334 while that of 
the odd is 9319 as before, but the state vectors assigned to these sectors are different since 
now we use (4.4) instead of (4.1) for the projection. In the unperturbed sine-Gordon model 
(A = 0) both projections can be done and the spectra compared to check the validity of the 
program, and we found a perfect match within the precision of floating point arithmetics. 

The spectrum does not change if we change the sign of y in either sectors, so the masses 
mf are even functions of A. This fact alone already precludes the existence of the linear 
correction (3.6), but we shall go through a detailed analysis in order to take into account 
the numerical uncertainties of the TCSA mass determination as well. We can also plot the 


9 




X 

mf 

m ! 

mf - nii (TCSA) 

mf — (3.6) 

0 

0.27230 

0.27231 

-0.00001 
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±0.010 

0.27154 
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0.00033 

±0.015 
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0.00002 
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±0.020 

0.26925 

0.26925 

0 

0.00066 

±0.025 

0.26750 

0.26746 

-0.00004 

0.00083 

±0.030 

0.26536 

0.26530 

0.00006 

0.00099 

±0.035 

0.26276 

0.26274 

0.00002 

0.00116 

±0.040 

0.25975 

0.25971 

0.00004 

0.00132 

±0.045 

0.25625 

0.25624 

0.00001 

0.00149 

±0.050 

0.25236 

0.25229 

0.00007 

0.00165 


Table 4.2: Masses measured from TCSA for (3 = ^4 , 6 — | in units of the unperturbed 
soliton mass M. TCSA mass measurements have an estimated error of ±0.00002. We also 
list the mass splittings measured from TCSA and compare it to the predictions (3.6). 


measured mass difference and the prediction from (3.6) as in figure 4.2, and the mismatch 
is obvious from the plot. 

We can also fit a function of the form 


mf = mi + a±x ± b±\ 2 ± c ± x 3 ± d±\ A 

to the measured masses in the x G [—0.05 ... 0.05] range. The fit results are 


a+ = 0 ± 0.00010 
b + = -7.606 ±0.018 
c+ = 0 ± 0.054 
d+ = -152.4 ±6.9 


a_ = 0 ± 0.000054 
b_ = -7.625 ±0.011 
c_ = 0 ± 0.033 
d_ = -149.5 ±4.4 


This must be compared to the prediction (3.6) which gives a± = ±0.0165367. We see that 
this value is two orders of magnitude larger than the uncertainty of the fit. Since the model 
in question is very deeply in the semiclassical value, even higher order quantum corrections 
cannot come to the rescue (we have seen in the previous section that they are at the level 
of a few percent). The conclusion is that the prediction (3.6) is not valid to the order (3 4 
which makes the conclusions of [7] based on these results unfounded (an analysis of the 
theoretical arguments in that paper is also given in the following section). 

We also analyzed the mass spectrum for [3 = This model has 8 breathers, so it 
is also very much in the semiclassical regime (/3 2 /8n = 0.10331...). Comparing FFPT 
with classical mass corrections again gives an estimate for quantum corrections at the few 
percent level. However, higher [3 means that the accuracy of TCSA is necessarily somewhat 
lower, but is still enough for our purposes. At truncation level 18, we found 7963 states in 
the even and 7978 states in the odd sector. The measured masses are summarized in table 
4.3, and the mass splitting is plotted in figure 4.3. 
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Figure 4.2: The mass splitting measured from TCSA compared to the semiclassical pre¬ 
diction (3.6) at (3 = ■yj? , 6 — |. The error bars indicate the estimated magnitude of 
truncation errors. 
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0.00137 

0.030 

0.35449 
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0.00005 

0.00165 

0.035 

0.35248 

0.35242 

0.00006 

0.00193 

0.040 

0.35011 

0.35011 

0 

0.00220 

0.045 

0.34747 

0.34751 

-0.00004 

0.00248 

0.050 

0.34443 

0.34445 

-0.00002 
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Table 4.3: Masses measured from TCSA for f3 = -yf , S — | in units of the unperturbed 
soliton mass M. TCSA mass measurements have an estimated error of order ±0.00005. 
We also list the mass splittings measured from TCSA and compare it to the predictions 
(3.6), 
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We can again perform a fit in the range x G [—0.05 ... 0.05] with the results 

a+ = 0 ± 0.00008 a_ = 0 ± 0.0002 
b + = -6.076 ±0.017 = -5.966 ±0.021 

c + = 0 ± 0.053 c_ = 0 ± 0.11 

d+ = - 63.3 ±6.5 d- = - 77.3 ±8.6 

The prediction (3.6) gives 

a± = ±0.0275751 

which is again ruled out by two orders of magnitude. 


5 Theoretical arguments 

In the previous section we demonstrated that TCSA data are strongly incompatible with 
the semiclassical analysis of the mass spectrum in [7], while they are in full accordance 
with the analysis of our earlier paper [4] based on form factor perturbation theory. 


5.1 Issues concerning the semiclassical form factor method 

Considering that the prediction (3.4) is only semiclassical, one may hope that higher order 
corrections can adjust it to be compatible to the numerical data. This hope is unfounded, 
however, since we have seen that quantum corrections can be estimated to be at most a 
few percent at the values of f3 we considered. Therefore we must conclude that (3.4) is not 
correct even to its lowest order. 

This is not surprising in the light of the fact that the leading order term (3.5) is of 
order /3 4 . The semiclassical computation presented in [7] is simply not valid to this order. 
Quantum corrections to soliton form factor and soliton mass should be taken into account 
before extracting the breather mass from the pole of the form factor. To order /3 2 , however, 
(3.5) is in accord with form factor perturbation theory that predicts no linear correction 
to the masses. 

We would also like to discuss briefly the mass corrections predicted in [7] for the case 
5 = 0. They obtain 


= ±- 


8y7^o 


— sin n^— + n log 


4). + 0(A «) 


Expanding in /3 


5rr^_ 

M 0 


± 



6144 



256 pi) 




(5.1) 


The classical limit 


6rr^_ 

M 0 
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is correctly reproduced by this formula. However, the f3 A term cannot be trusted because 
of the same reasons as above. 

Expansion of the FFPT mass correction (4.2) produces terms of all orders in (3 2 . It is 
obvious that a (3 2 term is also generated from (5.1) when we replace the classical parameters 
Ao and Mo by the PCFT coupling A and the quantum soliton mass M since the coupling 
and the mass are redefined at each order in the f3 2 expansion. For this reason alone the 
terms in (5.1) that are higher order in (3 cannot be compared directly to (4.2). 

However, there is a deeper problem with (5.1): the presence of the nonanalytic (loga¬ 
rithmic) dependence in Ao would suggest that although the masses have the correct limit 
when Ao —>► 0, the Ao ^ 0 case cannot be treated as a perturbation of the Ao = 0 model. 

In contrast, both conformal perturbation theory and form factor perturbation theory 
shows a rather different behaviour. Conformal perturbation theory produces a series ex¬ 
pansion in A and //, which contains no ultraviolet divergences at any order if (3 2 < 47t. 
There does not seem to be any way to produce nonanalytic behaviour in the coupling A in 
the regime of small (3\ even a partial resummation over /i would leave us with an analytic 
expansion in powers of A despite the fact that it may produce nonanalytic terms in // 
(due to infrared divergences present in the perturbation around the conformal field theory 
describing the massless ultraviolet fixed point). This partial resummation, on the other 
hand, is just expected to be identical to form factor perturbation theory (which is free 
of the infrared divergences since it expands around a massive field theory) which gives a 
first order correction in A with no log A terms, and matches very well with the numerics of 
TCSA as shown in subsection 3.1. On the other hand, since the semiclassical soliton form 
factor calculation in [7] is not self-consistent to order /l 4 , it cannot be used to argue for 
the presence of log A corrections. Before attempting any comparison to TCSA results, one 
must compute all the relevant loop corrections and suitably renormalize Ao, M 0 to A, M, 
which is out of the scope of the present work. 

We remark that even if all of the above is carried out, it could prove very difficult 
to provide any conclusive evidence against the existence of A log A corrections of order 
/ 3 A , since these cannot be separated easily from the already present A corrections, which 
are of order (3 °. However, the agreement of the measured mass splitting (illustrated in 
figure 4.1(b)) with the FFPT result (4.2) is very suggestive: the linear coefficient can be 
determined to a relative accuracy of 5 x 1CT 4 , and we still have not exhausted all resources 
to make the TCSA even more accurate: increasing the cut level is still a de fini te possibility 
(although it increases the required computer time very fast, with around the 3.2th power of 
the dimension of the Hilbert space), and extrapolating the cutoff dependence of measured 
masses to infinite cutoff is another one. 

5.2 FFPT and neutral states in the 5 = f case 

In the paper [7] the authors also produce a scenario in the 5 = | case to support the 
presence of the linear mass correction term (3.4). Let us recall the main line of their 
argument. 

It is known that the breathers arise as bound states in the soliton-antisoliton (ss) 
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channel, with alternating (7-parity. Namely, the nth breather B n with n odd/even arises as 
a bound state in a channel where the ss wave-function is odd/even under charge conjugation 
(interchanging soliton and antisoliton), respectively. 

The authors argue that such a situation is very unique to the sine-Gordon model and 
holds only if the non-integrable coupling A is not present. They claim that there exists 
another breather state of “wrong” C-parity and the n breather is actually a doublet B± 
(the upper index denoting (7-parity), but the states of wrong parity B~^~^ n decouple from 
the model which means that the corresponding three-particle couplings have the behaviour 

f$ = 0 for n odd /J 71 = 0 for n even 

so that only one pole appears in the ss channel for any n, with alternating (7-parity, since 
the other residue vanishes, and therefore the 2x2 ss scattering matrix degenerates to a 
one-dimensional projection at each pole. 

They continue by demonstrating that a modified form factor perturbation theory, ad¬ 
justed to take into account this degeneracy produces a linear correction to the breather 
mass of equal magnitude and opposite sign to the two breather copies, thus making the 
picture appear consistent with their semiclassical result (3.4). For h-folded models where 
the periodicity of the held is determined so that vacua are only identified after k periods 

2vr 

^ = ^ + k J] 

they predict 2k copies of each breather, twice as much as the h-folded model has, half of 
which is claimed to decouple from the model. 

There are, however, several problems with this scenario. Extensive investigation of 
the spectrum of sine-Gordon theory has been carried out using the TCSA approach. No 
sign of these extra states was ever detected, although the Hilbert space was guaranteed 
to be complete by taking all states of the ultraviolet conformal held theory into account 
(truncated by an upper energy cut, but all the experience with TCSA shows that this 
cannot be the reason for missing these states). Varying (3 down from the free fermion value 
V^br one can observe directly how the breather states appear. The result can be stated 
very simply in the case of the 1-folded model. We can identify the ss states in the spectrum 
by fitting their finite size energy dependence using Bethe Ansatz and the exact S'-matrix 
[13], and since the S'-matrix eigenvalue in the odd/even parity channel is different, we can 
even identify the parity of the states. When (3 crosses the threshold corresponding to the 
appearance of the n breather, a single ss state crosses below the 2 M limit (where M is 
the soliton mass), and this is a state of (7-parity (—l) n : new breather states appear one 
by one, and always with the right parity. There is no trace of the extra states, and this 
carries over to the Axfolded case as well: at each threshold value of (3 we always observe k 
copies of the appropriate breathers appearing in the spectrum (cf. [10] for some examples 
of spectra; we have produced many of them for different values of (3 , of which the published 
paper naturally contains only a few representative examples.) 

In the above argument, there is one possible loophole: TCSA uses periodic boundary 
conditions, and it may be argued that these exclude some states from the spectrum or 
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produce some identification between states, leaving only half of the breather states that 
are actually there in infinite volume 2 . However, even in this case one expects the mass 
spectrum extracted from TCSA to display the correct dependence on the couplings, and 
so the fact that our numerical results present such a strong case against any linear term 
(or, in fact, any odd power) in A remains a very compelling evidence against this scenario. 

The paper [7] also brings up an analogy with the theory RSOS 3 (the factorized scat¬ 
tering theory of the tricritical Ising model). The corresponding S'-matrix has no bound 
state poles, but such poles can be introduced by dressing with appropriate CDD factors. 
RSOS 3 has three vacua which can be illustrated by the potential 



The middle well of these is flanked by two others, and a simple calculation of the 
bootstrap spectrum reveals two breather bound states over it. These can be described as 
bound states of a kink going from the middle vacuum to one side and of a kink coming 
the way back. This seems to be analogous to the sine-Gordon vacuum structure, which is 
periodically infinite, but any chosen vacuum has again two neighbours. 

The analogy is flawed, however, because in sine-Gordon theory due to the periodic 
vacuum structure the boson is actually compactiked on a circle. In the standard (1-folded) 
case, the identification is 

2vr 



and therefore all vacua are identified. It does not make much sense to talk about two vacua 
flanking a middle one as in the RSOS 3 case where all three vacua are physically distinct and 
show up as 3 different energy levels in a finite volume (e.g. TCSA) spectrum, which become 
exponentially degenerate in the large volume limit. The 1-folded sine-Gordon model, on 
the other hand, has a unique ground state even in the large volume limit. In the fc-folded 
model there are k different vacua (with the topology of a circle), which is also manifest 
in the finite volume spectrum [10]. This vacuum structure is corroborated both by TCSA 
studies and by the exact description of the finite volume spectrum using the Destri-de Vega 
nonlinear integral equation (NLIE) [14] (see also [12, 13, 10]). It accounts readily for the 
k copies of each breather, but there is no additional doubling in the spectrum: the physics 
of these vacua is rather different from that of the RSOS 3 case. 

Our final conclusion is that the arguments given by the authors of [7] for these additional 
states are not well-founded, while the numerical data show that the dynamics of the model 


2 This issue was raised by G. Mussardo in a private discussion. There are known examples where the 
TCSA spectrum does not include all states of the infinite volume theory because of the selectivity of 
boundary conditions, however they are not very useful analogies of double sine-Gordon theory, similarly 
to the case of RSOS 3 discussed in the main text. 
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can be described very well by the usual form factor perturbation theory applied in [4] 
without any need for adjustments. 


6 Conclusions 

The numerical analysis has shown with quite large precision that the conclusions of [7] 
regarding the mass spectrum are untenable, more precisely the linear mass corrections at 
(5=1 predicted by them are ruled out by two orders of magnitude. The main issue seems 
to be that they trust the results from the semiclassical form factor approach to an order 
at which it is simply not consistent without performing the appropriate loop corrections 
as well. We also considered their theoretical arguments and have shown that they do not 
stand up to closer scrutiny. 

In contrast, the numerics is fully consistent with our original conclusions about the 
mass spectrum in [4], Obviously, this must not be taken as a complete verification of these 
results. It is always possible that there are effects that the numerics misses because of 
its finite precision and otherwise limited scope but which may be relevant for theoretical 
understanding. For the present we can only say that under closer examination the picture 
presented in [4] holds np. That picture is also consistent with many theoretical expectations 
and accumulated experience from the long study of two-dimensional quantum held theory 
(the ’folklore’), but certainly this is not a water-tight proof either. 

We conclude by noting that the double sine-Gordon model is an interesting quantum 
held theory in its own right, and we hope that continued work will shed more light on its 
workings. 
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